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GEOMETRICAL  OPTICS  OF  ANGULAR  STRATIFIED  MEDIA 

ABSTRACT 

Some  aspects  of  the  geometrical  optics  of  angular  stratified  media 
are  considered.  A  similarity  property  of  the  light  paths  enables  one  to 
draw  some  general  conclusions  about  the  behaviour  of  light  rays  traversing 
a  cone.  An  application  to  the  interferometric  method  of  observing  super¬ 
sonic  flows  shows  that  a  conical  flow  test  is  valid  including  the  refrac¬ 
tion  effect. 
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INTRODUCTION 


By  the  methods  of  geometrical  optics,  the  paths  of  rays  of  light 
can  be  found  by  integrating  a  system  of  ordinary  differential  equations. 

The  integration  problem  reduces  to  a  quadrature  in  the  two  special  cases 
when  the  isotropic  medium  is  (i)  plane  stratified  and  (ii)  spherically 
stratified.  In  these  cases  the  index  of  refraction,  n,  is  a  function 
of  a  single  variable,  in  (i)  a  Cartesian  coordinate  and  in  (ii)  the 
distance  from  a  point.  In  these  cases  it  is  possible  to  derive  a 
number  of  properties  of  the  rays  without  even  specifying  the  functional 
form  of  n,  [^l~fand  [2~j  .  In  particular,  the  light  rays  remain  in  a 
plane  so  that  these  two  cases  are  essentially  two-dimensional  problems. 

In  this  paper  we  consider  media  called  angular  stratified  media,  with 
the  property  that  n  is  a  function  of  the  polar  angle  with  respect  to  some 
line.  Here  the  two-  and  three-dimensional  problems  are  essentially  different 
however  they  both  have  a  "similarity  property".  This  similarity  property 
was  called  to  the  author's  attention  by  Dr.  J.  H.  Giese. 

In  neither  the  two-  or  three-dimensional  cases  was  it  found  that  the 
integration  could  be  reduced  to  a  quadrature.  But  some  general  conclusions 
regarding  the  light  paths  have  been  obtained.  In  the  two-dimensional  case 
the  problem  can  be  reduced  to  a  quadrature  plus  the  solution  of  a  first 
order  differential  equation. 

Consideration  of  angular  stratified  media  arises  naturally  in  the 
study  of  supersonic  conical  flow  e.g.  flow  over  a  cone  or  a  wedge  ( Taylor - 
Maceoll  or  Prandtl-Meyer  flow).  Since  optical  methods  of  observing  flows 
are  quite  common,  the  study  of  angular  stratified  media  is  not  entirely 
academic.  The  paper  concludes  with  an  application  to  the  interferometric 
method  of  observing  supersonic  flow  over  a  cone. 


THREE-DIMENSIONAL  CASE 


The  path  of  a  ray  of  light  in  a  medium  of  index  of  refraction  n 
can  be  written 


d(ndr/ds)  /ds  =  grad  n 


(1) 
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where  r:  (x,  y,  z )  is  the  position  vector  of  a  point  on  the  ray  and  s  is 
the  arc  length.  Equations  (l)  are  not  independent;  in  addition 

(dx/ds)2  =  1  (2) 

Suppose  n  is  a  homogeneous  function  of  degree  zero 

n(\r)  =  n(r) 

i.e.,  n  is  constant  on  radial  lines.  Then  grad  n  is  a  vector  whose 
components  are  homogeneous  functions  of  degree  -  1.  Then  it  is  easy  to 
see  that  (l)  and  (2)  are  invariant  under  the  similarity  transformation 

r— »  Xr 

s  — ►Xs  ( 5 ) 

This  is  the  similarity  property  of  the  rays.  An  angular  stratified 
medium  is  one  for  which  n  has  the  homogeneity  property  stated  above.  For 
example,  for  an  angular  stratified  medium  which  is  also  rotationally 
symmetric 

n  =  n(Vx+y/z)  (4) 

where  z  is  the  axis  of  symmetry. 

For  definiteness  consider  the  following  optical  problem.  A  cone  of 

half- angle  0,  vertex  at  the  origin,  and  axis  along  z, 

2  2  2  2 
x  +  y  =  z  tan  0 

is  imbedded  in  a  space  of  constant  index  of  refraction  N,  (See  Fig.  l). 
There  may  be  a  discontinuity  of  index  across  the  conical  surface.  Inside 
the  cone  the  index  is  of  the  form  (4). 
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A  parallel  beam  of  light  travelling  in  the  positive  x  direction  strikes  the 
cone.  The  ray  tracing  and  image  problems  are  to  be  considered.  In  the 
latter  the  mapping  of  points  in  some  image  plane  x  <0  onto  some  object 
plane  x  >0  is  to  be  studied. 

In  the  ray  tracing  problem  (l)  and  (2)  must  be  integrated  within  the 
conical  region  subject  to  initial  conditions  for  the  position  and  direction 
of  the  ray.  The  initial  direction  is  obtained  by  applying  Snell’s  law. 
Consider  the  points,  in  an  object  plane,  x  =  xQ-<0,  along  a  line  £ q  through 
the  origin,  slope  less  than  tan  0;  see  Figure  2. 


FIGURE  2 


OBJECT  PLANE 


The  rays  through  the  points  of  line  S( slope  tan  0)  strike  the  cone 
tangentially.  The  rays  through  the  points  of  £^  strike  the  cone 
along  a  generator.  Thus  the  initial  positions  of  these  rays  are  related 
by  the  similarity  transformation  r-*-X,r .  The  initial  direction,  from 
Snell’s  law,  is  invariant  under  this  transformation  since  this  direction 
is  determined  by  the  direction  before  refraction  and  the  normal  to  the 
cone  surface.  Therefore,  because  of  the  similarity  property  of  the  rays 
and  the  invariance  of  initial  direction,  the  light  rays  are  similar  curves. 
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If  and  are  the  position  vectors  of  two  rays  which  strike  the  cone 
along  a  generator ,  then  =  Xr^;  for  the  geometrical  path  length  at 
corresponding  points,  s2  =  Xs^.  Furthermore  the  optical  path  lengths 
andfiij  ,  where 

H  -J-- 

of  these  two  rays  are  related  by 

Kb x  [Li]  (5) 

Tracing  the  paths  of  all  rays  which  strike  along  a  generator  it  is  seen 
that  these  rays  lie  on  a  surface.  This  surface  can  be  generated  by  the 
motion  of  a  radial  line.  In  particular,  all  these  rays  terminate  along  a 
generator  and  therefore  all  leave  the  cone  in  the  same  direction.  We  have 

the  following  picture.  The  plane  sheet  of  rays  defining  in  the  object 

plane,  after  traversing  the  cone,  is  still  a  plane  at  some  angle  to  the 

original  plane.  Thus  in  an  image  plane  x  =  the  line  is  mapped 

into  a  straight  line  i.,  in  general,  not  through  the  origin.  Fig.  3- 


FIGURE  3 

IMAGE  PLANE 
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Under  certain  conditions,  e.g.,  n>N,  this  mapping  will  be  singular  over 
certain  regions  of  the  image  plane.  The  Jacobian  of  the  transformation 
will  be  zero  along  a  certain  curve,  the  envelope  of  the  lines  £ .  Between 
this  envelope  and  the  line  S  no  rays  are  imaged.  This  accounts  for  the 
"shadow  of  the  shock  wave"  in  a  shadowgraph  of  supersonic  flow  over  a 
cone.  Also,  for  n>N,  the  tangential  rays,  i.e.,  when  the  line  is  S, 
after  traversing  the  cone  do  not  emerge;  they  are  internally  reflected. 

This  will  also  be  true  of  a  small  bundle  of  rays  in  the  neighborhood  of  the 
tangential  rays.  The  last  two  statements  concerning  the  singularity  of 
the  mapping  and  internal  reflection  are  easily  proved  for  the  special  case 
n  =  constant  and  by  continuity  would  be  true  for  slowly  varying  n;  these 
results  have  been  verified  by  numerical  integration  of  (l)  and  (2)  for 
a  few  cases  with  variable  n. 

There  is  another  result  that  can  be  obtained  without  resort  to 
numerical  integration.  A  lower  bound  on  the  depth  of  penetration  of  a 
ray  into  the  conical  region  can  be  obtained.  In  QQ  it  is  shown  that, 
for  any  rotationally  symmetric  medium,  the  equations  of  the  light  rays  can 
be  transformed  into  the  equations  for  a  two-dimensional  problem  provided 
that  the  index  of  refraction  is  replaced  by  m, 

2  2  g  m  =[f2  " 

where  p  =  x  +  y  and  h  is  a  constant  for  a  given  ray.  This  constant 
can  be  written  as 

h  =  xq  -  yp 

where  p  and  q  are  the  optical  direction  cosines  with  respect  to  the  x 
and  y  axis.  For  our  parallel  bundle  of  rays,  evaluating  h  before  the 
ray  enters  the  cone 

h  =  -*yQ  N 

where  yQ  is  the  y  coordinate  at  the  piercing  point.  To  obtain  a  real 
value  of  m 

2  2  ,  2  2  2,  2 

p  >h  /n  >H  y0/nQ 

where  nQ  is  the  maximum  value  of  n  =  n(p/z).  For  example,  in  supersonic 
flow  over  a  cone  n  is  a  monotone  decreasing  function  of  p/z  so  that  n^ 
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is  the  index  evaluated  at  the  surface  of  the  cone.  In  the  few  numerical 
cases  considered  this  lower  hound  was  actually  found  to  he  quite  close 
to  the  point  of  deepest  penetration. 

TWO-DIMENSIONAL  CASE 


A  two-dimensional  problem  is  obtained  as  a  special  case  of  the  three- 
dimensional  problem  by  considering  a  parallel  sheet  of  rays  in  the  x  -  z 
plane.  In  this  case  introduce  polar  coordinate  r,  0.  Then  for  an  angular 
stratified  medium  n  =  n(0) .  It  is  evident  that  the  similarity  property 
holds  also  in  this  case.  Using  0  as  the  independent  variable  instead  of 
arc  length,,  the  path  of  a  light  ray  is  described  by  the  equation 


nr’/(r ' ^+r^) 


2.  2vl/2 


-  nr/(r,2+r2)1//2  =  0 


(6) 


where  '  =  d/d9.  Let  be  the  angle  between  a  radial  line  and  the  tangent 
to  the  path.  Then 

cos  =  r,/(r'2  +  r2)'1'^2 

sin  =  +  r/(r'2  +  r2)"*"^2 

Then  (6)  can  be  written 

-  (n'/n)cot  ^  =  +  1  (7) 

r'/r  =  cot  If  (8) 

The  integration  problem  reduces  to  solving  (7)  plus  the  quadrature  (8) . 

It  does  not  appear  possible  to  solve  (7)  without  specifying  n.  (The 

integration  problem  for  (7)  would  be  simple  for  special  functions  n, 
b0  \ 

e .  g . ,  n  =  ae  „ ) 


AN  APPLICATION 


In  the  interferometric  method  of  observing  supersonic  flow,  one 
tries  to  relate  the  density  in  the  flow  field  to  the  shift  in  a  fringe 
pattern.  Since  this  is  a  tedious  data-handling  problem,  it  is  useful 
to  devise  methods  of  getting  partial  information  from  a  photograph  of  the 
fringe  shifts.  One  such  method  is  a  test  for  conical  flow,  f~  3  H,  which 
is  a  result  of  showing  that  the  fringe  shift  is  a  homogeneous  function 
of  degree  one  in  y  and  z  In  H  this  was  shown  under  the  assumption 
that  the  light  rays  are  not  refracted  as  they  pass  through  the  conical 
flow  region.  The  same  result  can  be  obtained  allowing  refraction. 
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In  the  usual  experimental  arrangement  a  lens  is  placed  in  the  one 
beam  of  light  that  passes  through  the  disturbed  region.  The  lens  is 
placed  so  that  the  image  of  the  median  plane,  x  =  0,  is  photographed. 

Straight  lines  through  the  origin  of  an  object  plane,  are  mapped 
into  straignt  lines  through  the  origin  of  this  image  plane. 

The  fringe  shift  is  proportional  to  the  difference  in  optical  path 
lengths  of  two  interfering  rays,  one  of  which  passes  through  the  disturbance, 
the  other  does  not.  It  has  already  been  shown  that  for  a  ray  through  an 
angular  stratified  medium  the  optical  path  length  has  the  similarity 
property  (5).  It  remains  to  show  that  the  same  holds  for  the  optical  path 
length  of  the  interfering,  undisturbed  ray.  The  interfering  ray  can  be 
obtained  as  follows.  The  parallel  bundle  of  undeviated  rays  can  be  pro¬ 
jected  back  as  a  virtual  beam  (the  interferometer  actually  separates  this 
beam  from  the  deviated  rays).  Taking  account  of  the  lens,  it  is  necessary 
to  identify  an  undeviated  virtual  ray  with  a  deviated  one  that  appears 
to  come  from  the  same  point  in  the  median  plane,  see  Fig.  4,  which  is  a 
projection  of  the  three-dimensional  phenomena  onto  a  plane. 


y 


DEVIATED  RAY 


VIRTUAL  RAY 


MEDIAN  PLANE 


LENS 


FIGURE  4 
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The  appropriate  undeviated  ray  is  the  one  that  passes  through  Q.  Lay  off  a 
spherical  segment  with  center  Q,  and  radius  QP^ .  By  Fermat’s  principle 
the  optical  path  lengths  from  P^  and  P^’  to  the  point  at  which  inter¬ 
ference  takes  place  are  the  same  and  also  the  optical  path  lengths  from 
the  light  source  to  PQ  and  PQ*  are  the  same.  (Actually  the  interferometer 
is  usually  arranged  so  that  there  is  a  constant  difference  between  these 
optical  path  lengths  but  this  is  of  no  consequence  in  our  considerations.) 
The  difference  in  optical  path  lengths  between  P^P.,  and  P^P^  ’  must  be 
shown  to  have  the  similarity  property.  From  (5),  PQP1  has  this  property. 
Since  P^’  and  P  ’  are  determined  by  the  initial  and  terminal  points  of 
the  ray  P^P^ ,  the  coordinates  of  P^’  and  P^'  also  have  the  similarity 
property.  Since  the  virtual  beam  passes  through  only  a  region  of  constant 
index  N 

tP  ' P  ’  ~T  =  Ii  p  »p  » 

0  1  J  0  1 

lf  s  *  [Vi]  -  [w] 

then  5  has  the  similarity  property.  Finally  the  coordinates  in  the  median 

plane  (ym>  zm)  also  have  the  similarity  property.  Therefore, 

5 (Xy  ,  Xz  )  =  \5(y  ,z  )  (9) 

m  '  m  '  m  m'  '' 

This  is  the  basis  for  the  conical  flow  test  of  .  Note  that  it  is 

essential  in  proving  (9)  that  the  lens  be  used  to  focus  on  the  median 

plane.  For  any  other  plane  (9)  would  no  longer  be  true. 
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